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We consider the classical stochastic fluctuations of spacetime geometry 
induced by quantum fluctuations of massless non-conformal matter flelds in 
the Early Universe. To this end, we supplement the stress-energy tensor of 
these flelds with a stochastic part, which is computed along the lines of the 
Feynman- Vernon and Schwinger-Keldysh techniques; the Einstein equation is 
therefore upgraded to a so called Einstein-Langevin equation. We consider in 
some detail the conformal fluctuations of flat spacetime and the fluctuations 
of the scale factor in a simple cosmological model introduced by Hartle, which 
consists of a spatially flat isotropic cosmology driven by radiation and dust. 
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I. INTRODUCTION 



When a model of the Early Universe is developed, it becomes necessary to take into 
account the quantum nature of matter fields, even if for a wide range of energies, leading 
up to Planck's, the geometry of the Universe itself may be described in classical terms. 
The simplest and most studied way to couple quantum matter to a classical gravitational 
field is through the expectation value of the energy momentum tensor of the quantum fields 
IHJ^; however, the quantum fields must necessarily undergo fluctuations and it has 

been recently shown that these fluctuations are by no means negligible in many situations 
of interest [|,|| . 

To a certain extent, quantum fluctuations may be introduced in a classical model as 
uncertainty in the initial conditions. However, fluctuations play a subtler role when the 
semiclassical evolution, as it is in fact the rule, is dissipative [^p. In this case, the semi- 
classical gravitational field interacts in a nontrivial way with the environment provided by 
the quantum matter degrees of freedom P, leading, e. g., to anisotropy damping and loss 



of information regarding initial conditions |10-13|. The environment reacts back on the sys- 
tem in a way which is demanded by consistency of the underlying quantum theory. Under 
equilibrium conditions in static space times, the relationship between damping and back 
reaction becomes rather simple, and it is embodied in the so - called fluctuation - dissipation 
theorems Similar, if more complex, connections between damping and back reaction 

hold in nonequilibrium situations; as a rule, back reaction acquires a stochastic character 



in. 

In the particular case of semiclassical cosmology, we can be more precise about what is 
going on. Except in some very special cases, such as conformally flat universes coupled to 
conformal matter, or purely adiabatic evolution, the cosmological evolution leads to particle 



creation [16,17|. The phase relationships between the created particles are quickly lost; 
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even very weak interactions may accomplish this on very short time scales |]T8[. Under 
these conditions, particle creation is a dissipative process. The back reaction of the matter 
fields will then (loosely speaking) have a time reversible component associated to vacuum 
polarization phenomena (such as trace anomaly) and a time oriented component related 
to particle creation. The crucial point is that there is also a third, stochastic component, 
induced by the quantum fiuctuations of the matter fields. To keep the average effect while 
disregarding these fluctuations leads, as a rule, to an inconsistent theory. 

In recent years Langevin type equations have been proposed to describe the fluctuating 
back reaction of the matter flelds on the classical gravitational fleld, thus restoring the 
consistency of the semiclassical Einstein equations. The new semiclassical Einstein-Langevin 
equations predict classical stochastic fluctuations in the gravitational fleld 0,^,^,0. The 



derivation of these improved semiclassical equations is based, usually, in the Closed Time 
Path (CTP) effective action method and exploits the close connection that this effective 
action has with the Influence Functional introduced by Feynman and Vernon We 
should mention that a similar approach leading to Langevin type equations has been used 
in recent years to study the generation of primordial fluctuations in inflationary cosmologies 
2^^5| (Morikawa has independently arrived to similar conclusions regarding this subject 



||2q|), as well as the more general problem of stochasticity in effective fleld theory 

In this paper we study the stochastic fluctuations arising from coupling geometry to 
massless non-conformal flelds in two cases of interest. One of them is the study of the 
conformal fluctuations of flat spacetime. The second case is a spatially flat Friedmann- 
Robertson- Walker (FRW) spacetime with a classical source made of both radiation and 
dust. This model is simple enough that the fluctuations can be studied in some detail, 
yet it has signiflcant features of a realistic cosmology; it was proposed some time ago by 



Hartle [^, within the usual (Rosenfeld's) approach, to study the back-reaction effects to the 
creation of massless non-conformal particles in isotropic cosmologies. We shall not develop 



a full analysis, but rather focus on the new features brought by the explicit introduction of 
stochastic terms in the Einstein equations. 

The Hartle model has two perfect fluids, one of them with the equation of state of 
radiation and the other representing baryonic matter with the equation of state of dust. 
A dimensionless parameter, ^, measuring the relative amounts of baryons and radiation, 
which is constant in classical periods, is assumed to be always constant and with a value 
corresponding to the present universe (which is very small, ~ 10^^"^). A massless non- 
conformal quantum field is coupled to this system. The presence of the small portion of 
baryonic matter is essential to ensure particle creation, since in a radiative FRW universe 
the scalar curvature vanishes and the field cannot couple to the curvature. A perturbative 
expansion in terms of the parameter ^ is then seen to be equivalent to an expansion in a 
parameter that measures the deviation from the conformal coupling. Since the two degrees 
of freedom of the graviton field in a FRW background behave as massless minimally coupled 
fields 1^ Hartle's model provides a good testing ground for the study of the back reaction 
due to the creation of gravitons in cosmology. 

The model describing the conformal fluctuations of flat spacetime due to massless non- 
conformal fields can be easily obtained as a limit of Hartle's model. One just needs to 
impose that there is no classical matter source and that the cosmological scale factor is a 
small perturbation of unity. This model, which is considerably simpler than the cosmological 



model, may be of interest in connection with the semiclassical stability of fiat spacetime |^ 
(see also Ref. |^) under conformal fluctuations. 



We should stress that one of the main differences between the cosmological back- 
reaction problem considered in this paper and other cosmological back-reaction problems 
that have been considered in the literature is that the later usually deal with conformal 



fields in isotropic cosmologies [1 6,32, 10], or with conformal fields in cosmologies with small 



anisotropics p3| , P^pT| , p5| , pT| , p!^ or in cosmologies with small inhomogeneities [^,^. In the 



two last cases particle creation takes place because the anisotropies or inhomogeneities break 
the conformal invariance of the fields. 

Cosmological models where the back-reaction equations have been solved explicitly in- 



volve, usually, conformal fields in FRW cosmologies [^,^. In such models there is no 
particle creation and all quantum effects are a consequence of the fact that the vacuum ex- 
pectation value of the stress-energy tensor of the quantum field, subject to the trace anomaly, 
does not vanish. Generally the stress-energy tensor is also subject to a two parameter ambi- 
guity connected to the addition of two conserved terms quadratic in the spacetime curvature 
PJ40|1 . Such ambiguity cannot be resolved within the semiclassical theory and depending on 
these parameters one finds very different behaviors for the cosmological scale factor. Dif- 
ferent behaviors are also found as a consequence of the higher than two time derivatives 
of the scale factor that appear in the stress-energy tensor, this is a typical back-reaction 
phenomenon and the origin of runaway solutions. In this connection several authors have 
argued pi|j31| that solutions to the semiclassical back-reaction equations that deviate non 
perturbatively from the classical solutions should be discarded. 

When particle creation takes place, the back-reaction equations have been solved only 
numerically for the case of massive conformally coupled fields in FRW cosmologies In 
this case it is the mass that breaks the conformal symmetry of the problem. 

Here we are interested in the quantum effects due to the creation of massless non- 
conformal particles. Our dynamical back-reaction equations for the cosmological scale factor 
generalize that of Ref. in that a new stochastic source is included. Because we start with 



a OTP functional formalism our equations are real and causal and thus they admit an initial 
value formulation. Thus when our equations are averaged with respect to the stochastic 
source they should not be directly compared with the equations in Ref. [^], which were 



obtained using an in-out functional formalism. However the analysis carried out by Hartle 
on the behavior of the cosmological scale factor at large and very small cosmological times 



is still valid. For this reason we concentrate here mainly in the stochastic effects which are 
new. 

The plan of the paper is the following. In section II the actions for the massless non- 
conformal scalar field, the gravitational field, and the classical matter sources, which define 
Hartle's model, are described. In section III we give the stochastic back-reaction equations 
for the cosmological scale factor. The derivation of these equations is based in the CTP 
effective action describing the interaction of the classical scale factor with the quantum 
field, and the relation of this effective action with the Influence Functional of Feynman and 



Vernon ||2T|]. In section IV we consider the particular case of the conformal perturbations of 
fiat spacetime. In section V we consider the fluctuation in the cosmological scale factor in 
Hartle's model, and section VI summarizes our conclusions. 



II. HARTLE'S MODEL 



In this section we describe a model p8[ in which a massless non-conformal quantum scalar 
field is coupled to a spatially flat FRW cosmology with a classical source made of radiation 
and dust. Since we use a dimensional regularization technique some relevant terms in the 
classical action are given in n arbitrary dimensions. The cosmological model is described by 
the n-dimensional spatially fiat FRW metric. 



(ii.i) 



where 0(77) is the cosmological scale factor, rj is the conformal time and r/^i, = 
diag{—l, 1, ■ ■ ■ , 1). The classical action for the massless non-conformally coupled scalar field 
is 
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where R is the Ricci curvature scalar, 



+ is the parameter coupling the scalar field to the spacetime curvature; = -^zrj 
is the coupling parameter for a conformal field in n-dimensions. Since the metric ([1L1|) is 
conformally flat it is convenient to introduce a new rescaled matter field (j){x^) by (f){x^) = 

n — 2 

a^2- (^r])^{x^), and the scalar field action simplifies as, 



(f>{x^ 



(11.2) 



where now □ = 'q'^'^d^d,, is the fiat d'Alambertian operator. It is clear from this equation 
that if 1/ 7^ the conformally fiat symmetry is broken by the coupling of the scalar field. 

The action for the gravitational field, the Einstein-Hilbert action •S'^''^, needs to be cor- 
rected with a counterterm to cancel the divergencies that will come from the effective action. 
In our case it suffices to add a term quadratic in the Ricci scalar. Note that terms quadratic 
in Rfj^i, or the Weyl tensor could also be added but since these terms are not necessary for the 
renormalization we simply assume that their coefficients vanish. The divergent (in n = 4) 
counterterm is: 



2 n-4 



d^'x 



(11.3) 



where /ic is an arbitrary mass scale needed to give the correct dimensions to this action term. 

Since the only dynamical variable of the gravitational field is the cosmological scale factor 
a{r]) the gravitational action reduces in our case, when expanded in (ra — 4), to 

6V 



J dr] aa, 



(11.4) 
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(11.5) 



where = IGttGn [Gn is the Newtonian gravitational constant) and V is the volum integral 

To the above scalar and gravitational actions we need to add the action for the classical 
matter sources 5*^, which are a radiative and dust perfect fluids. This action is given by, 

S:!^[a] = -Vjdr] p,a, (11.6) 

where pb is a constant parameter related to the baryon energy density. It is connected to a 
similar parameter for the energy density of the radiation, through 



i= (11-7) 

Pr 

which is of order ^ ~ 10^^'' for the present universe. This parameter ^ measures the baryon 
to photon ratio. 

The action ( [11. 6| ) is justified because when it is varied with respect to the scale factor aijj) 
it reproduces the trace of the stress-energy tensor for a radiative perfect fluid with equation of 
state Pr = Pt-/3 (pr is the pressure and pr is the energy density of the radiation), and of a dust 
like baryon fluid of pressure = and energy density p^. When the dynamics of the scale 
factor is driven by the classical source only one finds that pb = PbO'~^ and pr = Pr^'^, then 
the relative amounts of baryons and radiation as defined by ipPb/pf^^ becomes constant and 
given by ( [II. 7| ). In fact, when the classical action ( [II. 6| ) is considered as the only dynamical 
source of gravity from the Einstein-Hilbert action ([II.4[) one gets 6a = {£p/2)pb or equivalently 
R = — (£p/2)T'^', where T'^' = —pb is the trace of the stress-energy tensor of a perfect fluid 
of dust (baryons) plus radiation. 

We note that when baryons are not present, i.e. when ^ = 0, the scalar curvature vanishes 
(i? = 0), and from the action ( [II. 2| ) for the scalar field 0(a;^) it is clear that this field behaves 
as a free field in flat spacetime. Therefore with the usual Minkowskian definition of a vacuum 
state (which corresponds to the conformal vacuum for $(a;^)) no particles can be created 



by the expansion of the Universe. Back-reaction effects are still possible due to the vacuum 
expectation value of the stress - energy tensor of the scalar field, which in this case depends on 
the two free parameters mentioned in the introduction. One of the parameters is associated 
to an action proportional to and the other to an action proportional to R^uR'^'^- The 
quadratic terms of the first parameter vanish when R = 0, but the quadratic terms of the 
second parameter are generally different from zero even when the scalar curvature vanishes. 
Thus all quantum effects in this case depend crucially on the second parameter being different 
from zero. We assume that such parameter is zero, thus = is a consistent solution to the 
semiclassical back-reaction problem. Since our interest is in the back reaction due to particle 
creation we will take ^ as a small non zero parameter, and we will consider the ua^R term in 



( p.2| ) as a perturbative term which gives a measure of the deviation from the radiative case 
R = 0. In practice we will compute the quantum corrections as perturbations in powers of 
the parameter z/, and this should be consistent with an expansion in ^. 

III. STOCHASTIC BACK-REACTION EQUATIONS 

In this section we compute the stochastic back-reaction equations for the cosmological 
scale factor a{ri) from the OTP effective action to one- loop order for the quantum field 



^(x^) in the FRW spacetime (|II.1| ). This action is evaluated using a perturbative expansion 



in powers of the parameter z/. The OTP effective action was introduced by Schwinger |43-45 



for its application in a curved background see Refs. ||T^j7|, ^ , p^ . The idea is to start with 
a generating functional from which one obtains expectation values, instead of the matrix 
elements one obtains using the generating functional of the ordinary in-out effective action 
P7| . The CTP effective action is the Legendre transform of such generating functional. The 
price one pays for this is that one needs to double the number of fields: the (+) and (— ) 



fields below. For the details of this paper we follow closely Ref. ||3^, see also Ref. 
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In our case we have two fields a{ri) and (j){x^) whose free actions are given by Sg[a; /ij + 
S^[a] from (|IT1]), ( [1T5|) and ( [1T6|) , and by = i / 0n0 from (|IT3), and an 



interaction action given by S'j„t[a, 0] = — | / d^x a^Rcj)^. The field a is classical, thus we only 
need to go up to one loop order for the field 0, which corresponds to the first order expansion 
of the generating functional in powers of h. We have 

+5„[a+,0+] - 5„[a",0_] - ^TrOnG), (III.l) 

where the fields (f)± correspond to the expectation values of and G, which is a ma- 
trix operator with 2x2 components, is the inverse of the classical kinetic operator 
A = diag{U-v{a+fR+,-{U-v{a-fR-)), obtained from ( [1X21) , with CTP boundary 
conditions (see |12[)- Since we are only interested in the effective action for the scale fac- 
tor a(?7), we set 0± = and expand the functional determinant in powers of the coupling 
constant; finally the renormalized CTP effective action reads: 

TcTP[a^] ^ S^Ja^] - Sfja'] + 5f^[a±], (III.2) 



where 



6 .2 9z/^ 



(III.3) 



and 



dr] r dr]' a(-) {r])N{r] - r]')A ( -) {r]'), (III.4) 
3 J-oo \aj \aj 



where we use the notation Af{ri) = f^irf) — f-ij]), {firf)} = f+irf) + f-{ji), for an arbitrary 
function f{ri). The kernels H(?] — 77'; /i) and N(77 — r]') are. 
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H(r/-V;/.) = -^/^e--"(-''')ln 
N(?7 — 77') = ':^~^iv ~ v')- 



p° + ie 



327r 



(III.5) 



Note that the parameter /i^ has been absorved into the new parameter Ji. 

The subscript IF in ^/^[a^] is used to indicate that this part of the CTP effective action 
is the influence functional action as defined by Feynman and Vernon PT|. That is, it gives 
the effect of the environment, the quantum field, on the system, the cosmological scale factor, 
which is the field of interest here. The identification between a part of the CTP effective 
action and the influence functional was suggested in Ref. M. 



A. Improved effective action 

The connection between the CTP effective action in the semiclassical context and the 
influence functional action introduced by Feynman and Vernon to describe the interaction 
between a system and an environment just mentioned gives an interesting new light on the 
semiclassical back-reaction problem. 

Let us now turn to the dynamical equations for the scalar field a{ri). In principle these 
equations can be derived from the CTP effective action as SrcTp/Sa'^\^±^^ = 0. A problem 
might arise from the fact that such effective action has an imaginary part, since Sj^p is 
complex, see ([1II.4|) . It should be clear however that the imaginary part of Sf^p will not 



contribute to the field equation derived in that form because Im^S^pla"^]) is quadratic in the 
difference of the (+) and (— ) fields and one finally takes = a~ = a. 

However, from the point of view of the system-environment, relation familiar in quantum 
statistical mechanics, the imaginary part of the influence action is related to the noise suffered 



by the system from the environment fluctuations This suggests an improvement to the 



semiclassical back-reaction equations that takes into account such fluctuations. This may be 
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achieved if we define the influence functional pT 



TiF[a^] = e'^^>['^ 1, (IIL6) 
and note that it may be written formally as 

^,^[a±] = |pj-p[^-]e^[«-(^f.[-*l)+6-/'^^^-WA(l)W], (III.7) 

with 

That is, performing the path integral in ( pi.7| ) with V [j] defined in (pi.8|) leads directly 
to ( |II1.6|) . If we interpret V [j] as a Gaussian probability distribution, the action in (p^lI.Tj ) 
may be seen formally as the action for a field a{r]) which is coupled to an external stochastic 
source j{r]). 

Thus, the influence functional ( [111.6| ) can be seen as the mean value with respect to the 
stochastic field j{ri) of an influence functional for an improved effective action Seff defined 
by, 

S,ff[a^;j] = S^Ja+] - S^Ja-] + Re{Sfp[a^]) + 6u J j{v)A (r^). (III.9) 



This comes from (|III.7|) and the addition of the gravitational and classical matter terms 



( [111.3|) . Now the field j{ri) will act as a stochastic source in the improved semiclassical 



equation when the functional derivation with respect to a'^{f]) is taken. This stochastic 
field is not dynamical, it is completely determined by the following relations, which may be 
derived from the characteristic functional, i.e. the functional Fourier transform of V [j], 

urn = 0, 
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Since the probability distribution is Gaussian the noise kernel is the two point correlation 
function of the stochastic field. In our case as one can see from (p.ll.5|) the noise is white. 



It is convenient for comparison with previous work p8| , p9| , to write the effective action 
( [111.9| ) in terms of dimensionless quantities such as x? Kx) ^iiid Cix) fo^' the conformal time 
T], cosmological scale factor a{ri) and stochastic source j{ri), respectively. These are defined 
by 

x-Itj". c-Itj^- 

We also introduce dimensionless frequencies u;, renormalization parameter /i and volume 
V instead of p", ]x and V, respectively, as: 



fil/2 , ^3/4 



= 6^/^£p/i, V = g^V. (III.12) 



Then the improved effective action ( [1I1.9| ) for the dimensionless scale field b{x) becomes, 



Z/2 



327r2 



+ :^C(X)A (X) }. (III.13) 
The stochastic source is defined, after ( [1I1.1CI| ), by (C(x)) = and 

The non-local operator k is defined by its action on a test function f{x) by, 

— e-'-^/zM/M, (III.15) 

-oo ZTT 

where f{u) is the Fourier transform of the test function, defined by /(tu) = f^oo^X^^'^^fix)^ 
and 
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h{uj) =^^n 



uo + is 



(III.16) 



This kernel is real and causal, as can be seen from its analytical structure. Note that 
this is in fact the main difference between the CTP approach and Hartle's approach. 

B. Stochastic back-reaction equations 

The dynamical equations for the cosmological scale factor in its dimensionless form. 



may be obtained by functional derivation of the effective action ( [III.13| ) as 



|T(5e//[&^;j]) 



b±=b 



0. After functional derivation and the identification of the 



and (— ) fields, the equation acquire an interesting form if we multiply it by h 

2 



- In 6 + 6— 

327r2 U; dx 



V b 



b d 

bdx 



(111.17) 



The effects of the quantum field are those proportional to the coupling parameter i/. The 
terms with In 6 come from the renormalization of the stress-energy tensor of the quantum 
field, see ( p.l.5| ), and there are non-local contributions in the terms with K[x',b/b]. The effects 
due to the creation of quantum particles, in particular the dissipation of the field b{x) are 
included in the non-local term. 

A point of interest is the role played by the parameter fi in the equation ( [ill. 17| ). This 
parameter is related to the local part of thus a change in this parameter is 

equivalent to a change in the coefficients in some local terms in the stress-energy momentum 
tensor of the quantum field. This parameter in fact is related to the two parameter ambiguity 
in the definition of such stress-energy tensor already mentioned. In our case one of the 
parameters has been fixed by our renormalization procedure and the ambiguity remains in 
the other parameter, see Ref. for more details. 
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When quantum effects are ignored, i.e. when we take z/ = 0, the back-reaction equation 
becomes 2h — ^ = 0. This is the classical equation when only the two classical ffuids are 
present; it admits the solution 

h = x + \^x', (in.i8) 

with appropriate initial conditions. When baryonic matter is not present, i.e. ^ = 0, the 
classical solution is b = x- If averaged with respect to C{x)y eq. (PII.17| ) also admits b = x 



as a solution, which is in agreement with Hartle's solution. As remarked in the introduction 
this is expected as in such case no particle creation takes place. Had we introduced local 
terms in the action of the type R^^R^^ , this would not be necessarily the case. 

IV. CONFORMAL FLUCTUATIONS IN FLAT SPACETIME 

The semiclassical back-reaction equation (P^ILlTp was derived for a non-conformal quan- 
tum field in a universe filled with classical radiation and baryonic matter, using perturbation 
theory around a conformal vacuum state. The classical contributions are in the first two 
terms of the equation i.e. 6^ — ^b. These two terms are obtained from the classical stress- 
energy tensor and different classical sources will lead to different terms. Thus by changing 
these first two terms in ( p.ll.l7| ), the equation may be used in different situations of interest. 
Besides the case of cosmological interest described in the previous section, a case of obvious 
interest is that of the conformal fluctuations of flat spacetime. 

To obtain the equations for the fluctuations around flat spacetime, one simple takes ^ = 
(no baryons) and notes that 6=1 (or a = 1) is the conformal factor for flat spacetime. Thus 
if we take 

Kx) = i + 9ix) (iv.i) 

in ( |111.17| ) with ^ = and consider the linear order in g{x) we find an equation for the 
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conformal fluctuations of flat spacetime, interacting with a massless non-conformally coupled 
scalar fleld. Equation ( |111.17| ) in this linear approximation simplifles considerably, it can be 
integrated twice and becomes 

,,2 
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i^lx; g] + ^Cix) = Aox + Bo 



(IV.2) 



where Ao and Bo are integration constants, Bo gives just a global rescaling of the conformal 
factor and Ao gives a linear expansion of the scale factor. This linear term is pure gauge, 
it represents an inflnitesimal coordinate change of t — > t + ^t^, where t is the Minkowski 
time. We can take Ao = Bo = a.s our background spacetime or, alternatively, deflne a new 
perturbation of the scale factor G{x) = 9{x) ~ ^oX ^ Bo, whose equation is 

,2 



(IV.3) 



The classical solution corresponds now to G{x) = (no perturbation). The quantum effects 
here reduce to the two terms which depend on u in ( [IV. 3| ), the flrst, of order u"^, is a non- 
local term which is linear in G{x)j and the second is the external stochastic source C(x)- The 
flrst includes the effects of particle creation, and the second drives the classical stochastic 
fluctuations of the conformal factor. To solve this equation we may take its Fourier transform, 

C(^) 



V 



where h(a;) is given in ( [1II.16|) . From eq. ([III.14| ) we know that 

V 



(IV.4) 



96 



and so 



{G{uj)G{uj')) 



?,MV 



1 + 
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h{uo) 



5{uj + u') , 



which leads to the real time correlation 
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<G(X)G(X')>^^/^ 



-iui{x-x') 



1 + 
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h{Lj) 



(IV.5) 



As 1/ — > 0, the conformal stochastic fluctuations become white noise, as it could be ex- 
pected. Indeed, we could Taylor expand inside the integration sign to develop the correlation 
function as a formal power series in v. An important problem, however, of this approxima- 
tion is that the fluctuation effect (in C(x)) and the dissipation effect (in k,\x', G]) are treated 
at different orders of approximation and this means that the fluctuation-dissipation relation 
cannot be fully appreciated. In fact, one would expect that as in the Brownian motion the 
scale factor driven by the stochastic term will loose energy by dissipation due to particle 
creation and that flnally a sort of equilibrium between the two effects will be reached. Here 
the stochastic effect dominates the behavior of b{x)- 

We may estimate the correlation time at flnite by rescaling uu = (87r/i/)s , we flnd 

r°° ds g-i87ri/-is(x-x') 



<G(x)G(x')) = ^£ 



27r 



1 2 



2 „4 



The integral is controlled by the small s range. Recalling that the function a;ln(x/e) has a 
minimum at a: = 1, we expand 



A 



where A = e{^T:/ ^vY, and we have used that A is a large number. In this approximation 



1 + s^n 



IIP 



+ 7rV~l + A(s2-^) +7rV 



A + TT^ 



/ - 28" + 1. 



The denominator has poles at 



A + 7r2 



±iVA + 7r2 - 1 + 1 



which converge towards 
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A + TT^ 



-1/4 ^ , X "^TT 

exp(l + 2A:)-, 



(A; = 0,l,2,3) asA— >oo. Choosing the contour of integration, we find 

z/2 e 



{G{x)G{x')) 



7687r\/ y/2 



cos 



eix-x'l- 



TT 



exp [-e|x-x'l 



where 



e 



4^2- 



TT 



(IV.6) 



We see that at finite u the correlation decays exponentially with a correlation time 0^^ ~ \/z/. 



V. FLUCTUATIONS IN A UNIVERSE WITH MATTER AND RADIATION 

We turn now to our main problem, which is the analysis of the fiuctuations of the cos- 
mological scale factor in Hartle's model, as described by (|111.17] ). 

The analytic solution of the stochastic integro-differential equation ( [111.17|) is not possible. 



but since ^, the parameter that gives the baryon to photon ratio is assumed to be very small, 
it makes sense to find solutions linearized around a radiative universe. Thus we follow Hartle 
28| and look for solutions of the type 



Substituting this into (pi.lT] ), taking only terms linear in ^, the equation can be integrated 
twice and becomes 

^ + TT^f-lnX-«:[x;-]U^- = ^ + A^x + Bo, (V.2) 
327r2x \X X J 2\/^x 4 

where Ao and Bo are integration constants which may be taken as Ao = Bo = 0, without 
loss of generality. The term gives the expansion corresponding to a matter dominated 
universe and comes from the second term in ( p.11.17] ). The terms with u are of quantum 



origin. Equation ( |V.2| ) is the dynamical equation for the perturbation g{x) of the scale 



factor around a radiative classical solution. 



This equation can now be directly compared with eq. (3.9) in Ref. |2^. In our case we 



have the external stochastic source ( which accounts for the fluctuations of the quantum 
stress-energy tensor and was not considered previously, but if we take the mean value with 
respect of the resulting equation is similar to Hartle's equation. The main difference here, 
as we pointed out earlier, is that the equation is real and causal. But the basic structure 
of the equation is similar and the conclusions that Hartle's draws in his analysis apply also 
here. 

In our case however we want to concentrate in the stochastic terms (which are new) and 
the dissipative terms, and how they contribute to the dynamics of g{x)- It is convenient to 
consider separately the behavior of fluctuations far away from and close to the cosmological 
singularity. 



To study the behavior of the solutions to the Einstein - Langevin equation, we first 
decompose g into the deterministic part g^ and the stochastic component gs- The behavior of 
the deterministic part follows Hartle's original analysis leading to g^ ~ x^/4+0(lnx/x^) 
for large x- 

As for the stochastic part, it is convenient to adopt as djTiamical variable / = gs/x^ 
rather than gs itself. The variable / satisfies 



A. Fluctuations far from the singularity 



dx XX (X - X ) f ix ) + {/ In X - [x\ /]} 



2Vi 



V 



c, 



(V.3) 



Instead of solving directly for /, it is convenient to look at the propagator 



G{x.xl) 



Sfix) 
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which obeys the equation 



f dx" XX" ix - X") Gix", X') + ^ {InxGix, x') - ^ [x; G{x, x')]} = -^Ax - x') , 



The main formal difference between this case and the conformal fluctuations of flat space 
time lies in that this equation is not time translation invariant. However, as in the previous 
case we expect the propagator to decay exponentially on the difference u = x ~ x'y with 
a weaker dependence on some 'center of mass' variable, which we take for simplicity to be 
X = {x + Our problem is then to disentangle the fast variable u from the slow one X 

p9| , |50[| . Now observe that, since we expect G{X, u) to be exponentially small unless u < X, 
we may approximate 



(V.4) 



lnxG(x,x') ~lnXG(X,M). 



On the other hand, let us write 



K [x; x')] = / rfx" Kx - x")Gix", X ) 



and set x" = x' + v io get 




Assuming that G depends weakly on X, we obtain 




Handling the classical term in the same way, 



p dx" Xx" ix - X") G{x", X) ~ X' f dv [u - v)G{X, v), 



we obtain the translation invariant equation 
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/ dv {u-v)G{X,v) 



327r2 

where X plays the role of a parameter. The solution reads 



{InX G{X,u) - k[u;G {X,u)]} = -7^S{u), 



a (X, u) = 



2ViJ 27r 



where 



G{X,u) 



X^ 



+ 



In 



e-^<^^G{X,uj), 



00 + ie 



(V.5) 



remarkably similar to its flat space counterpart. The self correlation is given in terms of the 
propagator as 

(/ (X) / (XO) = / dx" G(x, x")G{x': X"), 



which, with the same degree of accuracy, becomes 



(/ (x) / (xO) 



dio 



-ii^ix-x') 



\G{X,uj) 



We can repeat for this self correlation the same analysis than in the previous section, provided 
we replace {I'/X) for i> and ^X for /x, to get 



{9s (x) 9s (xO) 



cos 



©clx-x'l- 



TT 



exp[-ec|x-x'|]. 



(V.6) 



where ©c = XQ. As a result, the correlation conformal time 0^^ scales like ^JvjX. Of 
course, the approximations we have made assume that this correlation conformal time is 
small compared to X itself {X > i/^/^). In the region of validity, we see that the effect of 
expansion is to shorten the correlation time. 



B. Fluctuations near the singularity 



To obtain the form of the solutions valid at early times, we reason as follows. When x 
is small, the Fourier integral defining the nonlocal part of the Einstein - Langevin equation 
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is dominated by large frequencies a; ~ x ^. But since h (cu) depends only logarithmically on 
the frequency, we can then set h {u) = constant — h{x~^), to obtain 

« [x; / (x)] ~ - In [i^x] f ix) ■ 



The Einstein- Langevin equation then becomes 

,2 



X"9 + 



In 



X 



^X 



9 + 7^C 



167r2 \xoJ m 
where xo = Writing x = Xoe*^^^, we find 



X_ 
4 



(V.7) 
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+ d 



t / dt 



3t^ 



Defining a new variable v hy g = ie*^/^v, we get 



v" + -v' + 



32n\y _ 1 _ ^ 

1/2 4 



8^ Xo _ IStT^Xo 5tV4 



Again we decompose v into a deterministic and a stochastic part, Vd and Vg, respectively. 
The deterministic part admits a regular solution as t — > 



In the same approximation, we find 



t'2 

t 

t 



(more on ti later). Recalling that 



(V. 



we get, always in the same approximation, that 



{Vs (t) Vs (0) 



V 



/■*< ds 










t — 






L ^ 


t 
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where t< = min(t, t'). Note that if we assume that the stochastic source acts all the way 
from t = 0, we get an infinite amount of fluctuations. 

As a matter of fact, we do not expect our equations to hold up to the smallest times. First 
of all, our semiclassical universe presumably arises out of a fully quantum one at the quantum 
to classical transition. The semiclassical equation thus describes the model up to this time, 
but not earlier. Besides, we do not expect the form ([V.8|) of the noise correlation function 



to be true up to this time. The actual noise self correlation, which ought to be derived from 
quantum gravity, presumably has a finite correlation time (of the order of Planck's) and is 
not translation invariant, but depends on the distance from its arguments to the 'absolute 
zero of time' (to borrrow Misner's phrase ||51]]). Eq. (|V.8| ) is thus an approximation to be 
trusted only when the characteristic times are much larger than Planck's. The lower cutoff 
ti can be thought of as the earliest time when the model becomes meaningful. 

While we cannot make a definite prediction for ti, we can still observe that if it is 
small enough, then fluctuations build up fast enough that they may actually dominate the 
deterministic part over an stretch of time, namely, when 

,1,1 ^ 

Int < —= T7^- 

It is interesting to look at the behavior of a quantity of direct physical significance, such 
as the Hubble parameter Ti = a/a"^, rather than the metric itself. In terms of dimensionless 
quantitities (cfr. eq. ( [1II.11|) ) H = H/\/6lp, where 

H=--. 

dx 

Expanding b{x) = X + ^9{x) we can write H in terms of t { x = Xo^*^^^) and v as 



X 



1 + le-^/^ + (^^^ 
Xo \dt \ 2t ^ 



V 



Now, decomposing as above H = + Hs, vte obtain 
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which leads to 



X 



1 + 



;21nt + 1) 

















" 1 


1" 


\ vVx^ ) 3847rXo 







We see that, provided ti is small enough, the fluctuations may have an strong effect on the 
Hubble parameter. 

VI. CONCLUSIONS 



In this paper we have considered the stochastic, classical metric fluctuations induced 
by quantum fluctuations of matter flelds in two cases of interest, namely the conformal 
fluctuations of flat space time and a simple cosmological model, the former being mostly of 
interest as preparation for the latter. The quantum fluctuations act on geometry through the 
energy momentum tensor, which has a deterministic part, associated to vacuum polarization 
and particle creation, and also a fluctuating part, related to the impredictable aspects of 
quantum behavior; the introduction of the later leads to the formulation of the theory in 
terms of a so-called Einstein-Langevin equation. 

While the physical rationale for including such terms is clear (and earlier estimates point 
to their quantitative relevance PJ^) and the technical means to introduce them consistently 
into Einstein equations are by now well established ||8| j2^JT5|JT9|j20[| this is, to our knowledge, 
the flrst time actual non trivial solutions to the Einstein-Langevin equation are obtained. 
The reason why the solution to particular models lagged so much behind the formulation of 
the basic principles lies of course in the staggering complexity of Einstein theory, involving 
issues of gauge invariance, renormalization and nonlinear dynamics. The model proposed by 
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Hartle in Ref. is an ideal testing ground because it retains most of the physical features 
of a realistic cosmological model, while being simple enough to be tractable. 

Since in this last analysis massless non-conformal fields in FRW backgrounds are of inter- 
est mostly as a toy model for gravitons (for the application of similar ideas to inflationary 
cosmology, see the effects we have discussed here are at bottom quantum gravi- 

tational, and they are most relevant in the earliest stages of cosmic evolution. This adds a 
degree of uncertainty in the theory, associated to the present lack of understanding of the 
details of the quantum to classical transition; the best we can do is to localize our igno- 
rance in a few undetermined parameters, such as the time at which the initial conditions for 
semiclassical evolution are set and the value of the several renormalized parameters in the 
model (in particular, we have applied Occam's razor in setting to zero any parameter not 
specifically demanded by renormalization). 

Given these limitations, the main conclusions of our effort are that the Einstein-Langevin 
equation may indeed be solved, that a consistent picture of semiclassical evolution emerges 
and, most important, that this picture is significatively different from what earlier (non 
stochastic) semiclassical models have led us to believe (consider for example the ratio of the 
fluctuating to the deterministic parts of the Hubble parameter, as discussed in last section) . 
Clearly more work shall be needed before questions of direct cosmological impact (such as 
whether stochastic fluctuations may be instrumental in creating the homogeneous patches 
where inflation becomes a possibility, or rather work against the stability of those patches) 
may be addressed. We continue our research on this, in our opinion, most relevant issue of 
pre-inflation cosmology. 
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